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ABSTRACT 

za.sc* 

Considering the numerous and often differing theories available for deter- 
mining the behavior of compressible turbulent boundary layers a choice of pre- 
diction method is often determined by the available experimental data. This 
paper considers this problem in the light of the results of recent experiments 
in a number of different facilities and from flight. The data for skin fric- 
tion and heat transfer are generally from sharp flat plate and pointed cone 
models and are now available for very cold walls up to Mach numbers of approx- 
imately 9 - Mach number and wall-temperature effects from these experiments are 
compared to the more recent prediction methods. In addition, an evaluation of 
the maximum heat transfer and shear stress to be encountered on smooth flat 
plates and cones is made. Since boundary-layer s tructu re is an important con- 
sideration, boundary-layer v elocit y and temperature profiles obtained on hollow 
cylinders and nozzle walls are considered. 


INTRODUCTION 


A large number of methods for predicting the behavior of compressible tur- 
bulent boundary layers have accumulated over the years. Often the results from 
these methods differ widely and recourse to experiments must be made. In turn 
the methods themselves are often based on experiments so that an alternating 
procedure develops. In our case, we will bring to bear on this problem the 
results of recent experiments in a number of different facilities and from var- 
ious organizations. 5 


The results for heat transfer and skin friction that are presented here are 
in general restricted to the flat plate with a "sharp" leading edge and the 
pointed cone; however, boundary-layer profiles obtained on hollow cylinders 
and nozzle walls are also presented. Evaluating the data obtained on the flat 
plate and cone models is complicated by problems such as establishing the 
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virtual origin of the boundary layer, the Reynolds analogy factor, and the 
relationship of results obtained on cones to those obtained on flat plates; 
these problems are discussed in the presentation. 

The results allow a simultaneous assessment of wall-temperature effect and 
Mach number effect (Mach numbers up to about 9) in relation to the more recent 
prediction methods. Additionally, consideration is given to evaluating the 
maximum heat transfer and skin friction that is encountered on smooth flat 
plates and cones from very low to hypersonic speeds. 


SYMBOLS 


a .w. 

c F 

C f 

c f,i,x 

d 

f.p. 

h 

L(K) 

L(P) 

M 

n 

N Pr 

%t 


adiabatic wall 

average skin-friction coefficient 
local skin-friction coefficient 

incompressible value of turbulent Cf at value of where 

peak in Cf occurs (f^p) 

leading-edge thickness 

flat plate 

enthalpy 

Karman mixing length 
Prandtl mixing length 
Mach number 

inverse of the exponent in the power law for skin friction 
Prandtl number 

heat-transfer coefficient (Stanton number) 


^St,c heat-transfer coefficient with boundary-layer turbulent from cone 

apex 

Nst,i,c incompressible value of Stanton number at R v / 2 

%t,i,x incompressible value of Stanton number at value of R where 

peak in N gt occurs 

r radius of axisymmetric nozzle 
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R Reynolds number 

Rq Reynolds number based on momentum thickness 

R v Reynolds number based on local conditions and distance from peak 

shear or peak heating 

R, value of R 7 where peak shear or peak heating occurs 

t j p L 

Rj Reynolds number based on local conditions and distance from 

leading edge 

t time 

T temperature 

u local velocity in boundary layer 

U velocity outside boundary layer 

y distance measured from surface, normal to surface 

a angle of attack 

5* displacement thickness of boundary layer 

9 momentum thickness of boundary layer 

0 C semiapex angle of cone 

od exponent in power law for viscosity 

Subscripts: 

aw adiabatic wall 

c cone 

fp flat plate 

i incompressible value 

l local value 

p peak value 

t total value 

tc truncated cone 
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w 

wall 

00 

free stream 


DATA ANALYSIS 


Of importance in evaluating sets of data which cover a large range of con- 
ditions is the determination of a virtual origin, which is hopefully consis- 
tent. Among many assumptions one can use (a) the physical leading edge; 

(b) the start of transition [l] ; (c) the end of transition (say, the place of 
peak shear stress or peak heating) H H [“] ; or (d) match the momentum loss 
across an abrupt transition point and theoretically determine a hypothetical 
start for the turbulent boundary layer £ 5 ]] [jQ* 

The first three assumptions were checked against a large body of heat- 
transfer data obtained at the Langley laboratory and certain low speed experi- 
ments. The only one of these various assumptions of effective boundary-layer 
origin, which was reasonably consistent for the bulk of the data was assump- 
tion (c) which utilizes the peak shear stress or peak heating as the origin, 
thus agreeing with Coles early work Qf] . This assumption of virtual origin was 
used for the work reported in this paper. As a consequence of using this vir- 
tual origin, only data where a peak can be identified, or where transition 
close to a trip device occurs, can be used. Assumption (d) or momentum 
matching is most attractive from the standpoint of being the least arbitrary 
but requires a priori selection of a compressible theory for its utilization 
in general. 

Low speed heat-transfer data [jiT] are shown in figure 1 with Reynolds num- 
bers based on distance from the leading edge and from the assumed virtual 
origin (peak) . In the turbulent case the Karman-Schoenherr equation for local 
skin-friction coefficient modified by the analogy between heat transfer and 
momentum according to von Karman is shown for comparison purposes. (See 
appendix A for the incompressible formulas used.) The data generally show good 
agreement with the theory when measured from the peak location. Exceptions are 
the data at the lowest Reynolds numbers (i.e., closest to the peak) and the data 
where natural transition occurred. This result at the lowest Reynolds numbers 
is expected considering the shape of the experimental curve. However, the dif- 
ficulty with the natural transition data is believed due to the use of the 
indicated peak in the data which form a long shallow curve in the region of the 
peak. A relatively small forward correction to the peak location would put 
these data in good agreement with the other data and the theory. 

It might be well to point out here that the Reynolds analogy factor chosen 
to modify the skin-friction equation can have a significant influence on the 
level of the experimental data in the form of the %t/%t i ra ' tio ’ The 

Karman analogy factor used |jfj ^with Np r = 0-725, for this paper) is within 
±2 to ±3 percent of that given by the more recent papers of Deissler-Loeffler 
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(M = 0) J~£TJ and Tetervin QfJ tut is 9 to 11 percent lower than the value given 
by the Colburn analogy (%> r “2/3 with Np r = O. 725 ) for 10^ < R v < 10 7 and 
there are still greater differences at higher Reynolds numbers. 

Examples of some high-speed unmodified data are shown in figure 2(a). 

These data were obtained by Cary in the Langley 20-inch hypersonic tunnel on a 
flat plate 16 inches long and 11 inches in span. Data at three pressure levels 
and two angles of attack are shown. The highest angle of attack reduces the 
free- stream Mach number of 6 to a local Mach number of 4. The known effect of 
leading-edge thickness and unit Reynolds number on boundary-layer transition [lo] 
seem sufficient to explain the movement of transition with changes in pressure 
level and here one may note that station A and station B designated on the fig- 
ure are located one-half inch on either side of the model center line or only 
one inch apart. However, we are not so much concerned with transition per se 
as in simply knowning where the peak occurs. 

These same data plus data obtained at other angles of attack are pre- 
sented in figure 2(b) with Reynolds number based on distance from the peak 
heating value. In this case the experimental points have been divided by the 
values given by the modified Karman-Schoenherr equation at the corresponding 
values of R v . One can see that for R v greater than about 10^ the normalized 
heat-transfer coefficient is virtually independent of Reynolds number. The 
same result is indicated by data shown in figure 3 obtained in the Langley 
18 -inch variable density wind tunnel on the same model tested at a free- stream 
Mach number of 8 by Weinstein. 

Both skin friction and heat transfer have been obtained on a flat plate by 
Neal in the Langley 11-inch hypersonic tunnel and are shown in figure 4. In this 
case the model was 10 inches in span and 2k inches long; however, the tunnel 
test core width was 5 to 6 inches. Data is shown to a distance of 19-75 inches 
from the leading edge along the model center line. Visual studies using sur- 
face oil flow indicated some side edge effect near the trailing edge. Skin 
friction was obtained by the floating element technique and the heat transfer 
by using thermocouples in a thin skin. Different leading edges were used for 
the skin friction and the heat-transfer measurements, and even though both 
measured about 0.001 inch thick, undetected differences between the two leading 
edges may explain the change in transition Reynolds number between the two sets 
of data. This however does not explain the insensitivity of transition to 
pressure level. The heating data as a function of distance from the peak 
heating (R v ) shows much the same behavior as the experiments presented in fig- 
ures 1 to 3 although there is a tendency for this data to show a slight decrease 
for the highest values of R v . This may be attributable to the side edge 

effects previously mentioned. 

Thus far the data considered has been from flat-plate models; however, 
much of the data available is from cones. According to Van Driest [ll] the 
local skin-friction coefficient on a cone is the same as that on a flat plate 
at half the local Reynolds number on the cone. This rule however applies only 
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to cones where the flow is turbulent from the apex; in the actual case transi- 
tion occurs at a significant distance downstream of the apex and the turbulent 
boundary layer would be expected to behave more as if it were growing on a trun- 
cated cone. A transformation has been applied to the data obtained on cones to 
correct the results to the values that would be obtained were the flow turbulent 
from the apex of the cone. The details of this transformation are given in 
appendix B. 

Two sets of data obtained on cones in supersonic wind-tunnel flow [ 12 ] [l5j 
are shown in figure 5- Here only the data downstream of peak heating are 
shown. The heat-transfer coefficient has been normalized with respect to the 
incompressible heat-transfer coefficient, as for the data in figures 2 to 4, 
except that the incompressible value at half the Reynolds number (R v ) of inter- 
est is used, designated in this case Ng^ ^ c - For this and the remainder of 
the cone results, the data is reduced in terms of local values on the cone sur- 
face M- The amount that the original data (open symbols) lies below the 
modified data (filled symbols) indicates the reduction in heat-transfer coef- 
ficient that theory predicts to be due to the location of turbulent flow down- 
stream of the cone apex. The effect is seen to be significant. 

The results of three separate flight experiments Jl5] [l 6 j [if] are given 

in figure 6 . One notes that the original normalized heating coefficient data 
values (fig. 6 (a)) are generally increasing with increasing Reynolds number 
(R v ). This sort of effect is predicted by the theory taking the "turbulent 
cone truncation" into account. (One can note a general tendency of this type 
in the overall original data of figure 5 , but the restricted range of Ry and 

data scatter did not allow definite conclusions to be drawn.) Upon examining 
the individual series of datum points corrected to turbulent flow from the apex 
(fig. 6 (b)), the effect of Reynolds number on the Stanton number ratio is, at 
least, found to be significantly reduced. (Datum points close to the place 
where peak heating occurs are not included in the corrected data; these may 
readily be identified by comparing the two parts of figure 6 in the lower 
range of R v • ) 

Of considerable value is the hitherto unpublished data, from cone models, 
contained in figure 7* These data cover a considerable range of wall tempera- 
ture ratio and most importantly extend down to very low wall temperature ratios. 
The upper set in figure 7( a ) is derived from data supplied by the Republic 
Aviation Corporation from tests in their 36 -inch-diameter pebble bed heated, 
closed jet, hypersonic wind tunnel. ^ The lower set is an adaption of data sup- 
plied by the McDonnell Aircraft Corporation from tests that were run in the 
Cornell Aeronautical Laboratory 48-inch hypersonic shock tunnel. 5 The manner 
in which wall -temperature ratio was varied was different for these two sets of 
data. For the Republic data the supply temperature was kept roughly constant 
and the wall temperature was varied by cooling or heating the model. In the 
case of the McDonnell data, the wall temperature was always about room temper- 
ature and the wall -temperature ratio was varied by using a number of different 

^U.S. Air Force contract AF 33(657) - 9697 . 

5u.S. Air Force contract AF 04(694) -390. 
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supply enthalpies. Again the correction for "turbulent cone truncation" 
reduces the effect of Reynolds number on the Stanton number ratio (fig. 70 3 ))? 
though the effect of Reynolds number in the case of the Republic data is more 
than would be expected, it is not considered excessive. 


EFFECT OF WALL-TEMEERATURE RATIO ON HEAT TRANSFER 


In the previous section the data was considered from the standpoint of 
its proper reduction and interpretation rather than examining trends or evalu- 
ating predictions. Here we will examine the data in relation to the prediction 
methods. One of the important points to be resolved in the turbulent boundary- 
layer problem is the effect of wall-temperature ratio on the heat-transfer (or 
skin-friction) ratio. To investigate this problem in some detail, figures 8, 

9, and 10 have been prepared to cover the supersonic-hypersonic range and to 
compare flight and wind-tunnel results where possible. 

Data obtained in the supersonic range between Mach numbers of 3 and 4 are 
shown in figure 8. These data are average values of the data contained in 
figures 2(b), 5> 6(b), and [l8j. measureme nts contained in &<Q were 

taken on the cylindrical portion of a cone cylinder. The actual Mach numbers 
shown are 3*3 and 3*8. Where the Mach number of the test was different from 
these Mach numbers, the predicted incremental effect of Mach number on the 
heating coefficient ratio was added to the datum value. In this range of Mach 
numbers there is an insignificant difference between the various theories in 
the predicted increment of skin friction or heat transfer due to a small change 
in Mach number. The theories shown are the much used T-prime or reference tem- 
perature method using Monaghan f s constants as given in [19] and the more recent 
semiempirical formulation for skin friction of Spalding-Chi |pcTJ modified to 
heat transfer as given in appendix A. 


The trend though not the magnitude of the M = 3-3 data is predicted 
quite well by the T-prime method and the trend of the M = 3-8 data appears to 
be intermediate between the T-prime and the modified Spalding-Chi predictions. 8 
At both Mach numbers there is generally good agreement between the flight and 
wind-tunnel results. It may be noted that the T-prime method was computed for 
a nominal flight temperature at both Mach numbers. Computing the T-prime result 
for the lowest wind-tunnel total temperature (500° R) raises the prediction 
about 5 percent over the curves shown. As given, the Spalding-Chi method does 
not include an effect of total temperature level. The M = 3*3 data are con- 
siderably and consistently above both prediction methods - an unexplained 
effect. However, data obtained on wing panels during flights of the X-15 by 
Banner, Kuhl, and Quinn of the NASA Flight Research Center at a Mach number 
of 3 do not agree with the high level indicated by the M = 3*3 data in fig- 
ure 8. In the case of the X-15 data at a Mach number of 3; which is shown on 


^Coles results 00 for local skin friction at M = 3-7 with an adiabatic 
wall ^T w /T t « 0 . 9 ) give a value for CfjCf of about 0.54 thus fitting in 
quite well with the heat-transfer data. 
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the left-hand side of figure 9 "by "the shaded rectangle, there is good agreement 
with the modified Spalding-Chi prediction. As shown to the right in this fig- 
ure, X-15 data at a Mach number of 5 still agree best with the modified 
Spalding-Chi prediction, but the mean of the data appears to lie somewhat below 
the prediction. The value of %t,i use< ^ to normalize the experimental heat- 

transfer coefficients is from the modified Karman-Schoenherr equation as for 
the other data presented in this paper. 


Flight data are desirable but because such data generally cannot be 
obtained under as well controlled conditions as are possible in wind tunnels, 
doubts can always remain where such experiments are our major source of infor- 
mation. Thus much of the available experimental results from wind tunnels have 
been collected in figure 10 to aid in evaluating the wall-temperature effect. 
The data are in the low hypersonic speed range (roughly, Mach numbers from 5 
to 9 ) where the large temperature potential driving heat transfer aids in 
obtaining more accurate data than are generally obtained in the supersonic 
speed range in wind-tunnel flows. The results shown are average values from 

figures 2 (b), 3 , 4, and 7 (b) plus additional data from |l 8 j • M- 

and [ 23 ]. 


For the data shown in figure 10 which were obtained by the conventional 
thin wall thermocouple or thin-film techniques (solid symbols), good agreement 
is found with the trend of the Spalding-Chi predictions. The T-prime method 
significantly overpredicts the heat transfer at the lowest wall-temperature 
ratios. This overprediction is not apparent unless the wall-temperature ratio 
is considerably below 0-5 at M 85 5 and less than 0.3 to 0.4 at M 8 . The 


NOL and APL data |~2lj , jj 22 j, and QjjjJ (open symbols) generally do not agree 
with the other data. The Winkler-Cha formulation was designed to fit the 


NOL data, which it does quite well, but except for Hill's data considerably 
underpredicts the other data including the X-15 results previously presented 
(fig. 9 ) . The NOL heat-transfer results were obtained from inner wall- 
temperature measurements and the data of Hill were obtained from boundary-layer 
profiles; Winkler and Cha's results were on a flat plate, Lobb, Winkler, and 
Persh obtained their results in a wedge nozzle, and Hill's results were 
obtained in a conical nozzle utilizing nitrogen gas. 


Going back to the data represented by the solid symbols in figure 10, one 
notes that, although the modified Spalding-Chi method best represents the data 
of all the theories shown, it does somewhat underpredict the data. The under- 
prediction is generally from 5 to 15 percent with one point at Mach 8 about 
25 percent underpredicted. This discrepancy appears to be a function of, and 
increasing with, Mach number. It is perhaps best seen in the original presenta- 
tion of the data in figures 2(b), 3 , 4, and J. It is clear that we cannot at 
this point say that this represents a basic effect, for the problem of the 
proper form of the Reynolds analogy factor relating skin friction and heat 
transfer has not been solved. For instance, using Colburn's form of Reynolds 
analogy (discussed earlier) would bring the prediction into somewhat better 
agreement with the higher Mach number results but would result in poorer agree- 
ment for the lower Mach number data (fig. 1 results and fig. 2(b) M = 4) . 

There are possible, however, other explanations for this type of behavior. 
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Various turbulent theories are quite plastic in their behavior depending 
on the assumptions which go into their makeup. The predictions of average and 
local skin-friction ratios on an insulated flat plate as a function of Mach 

number by various theories [20], [24], [25] are shown in figure 11. Previously 

in this paper we have avoided such a presentation of theory because of the dif- 
ficulty of comparing data and theory on such a plot when, in addition to the 
conditions shown, the problem of the large range of wall-temperature ratios and 
stream conditions has to be added. Often in the past such plots on which are 
placed all available theories have been used as a horrible example to show the 
large discrepancies between theories. Here our intention is more to show the 
wide variations in prediction possible within a given theory. 

Generally the more complete theories show an effect of stream temperature 
level. In this case, this is illustrated by showing calculations for two limits 
of the exponent in the power law for viscosity, 00 =0-5^ which corresponds to 
a high Mach number - high temperature flow, and (JO = 1, an arbitrary upper 
limit which might correspond to a very cold flow. (For the tests presented in 
this paper cd generally varied between 0.8 and 1.) In addition, the analytic 
result depends upon the law assumed for the mixing length; either the Prandtl 
mixing length (L(P)) or the Karman mixing length (L(K)) are generally used. 

Li and Nagamatsu [25] also introduce a proportionality constant or "compressi- 
bility mixing length parameter" which determines the importance of the contri- 
bution to the shearing stress from density fluctuations. When the proportion- 
ality constant in the Li-Nagamatsu theory is zero, the result is the same as 

that given by the Van Driest theory . Since the exponent to is generally 
known, the most important determinations are that of the proper mixing length 
and of the contribution of the density fluctuations. The Spalding -Chi method 
is semiempirical and does not, as given, allow a choice of mixing length or 
values of to. The T-prime method, included for reference, does not have mixing 
length considerations. Figure 11 was prepared for a Reynolds number of 10?. 
Changing the Reynolds number by an order of magnitude, say, would change the 
curves shown but this would probably be a second-order effect. 

In the Li-Nagamatsu paper local skin friction and heat transfer and the 
effect of wall-temperature ratio are not treated directly; however, because of 
its similarity with the Van Driest analysis together with the increase in skin 
friction due to the inclusion of density fluctuation terms, it is possible that 
an extension to this theory may give a good prediction of heating at the higher 
Mach numbers if the Karman mixing length is assumed to be the proper one. This 
result can only be considered tentative and definitive answers will probably 
have to await experiments at still higher Mach numbers than have been presented 
here. If this result is correct, however, higher heat-transfer rates than are 
predicted by the Spalding-Chi method will have to be anticipated at high Mach 
numbers. 

The effect of wall-temperature ratio was one of our major considerations 
in the data in figures 8 to 10. How the local skin-friction ratios presented 
in figure 11 for the insulated wall case are affected by temperature ratio is 
shown in figure 12. Again the limits for the viscosity-temperature relation 
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exponent are taken as 0-5 and 1 (except for the Spalding-Chi method) and where 
applicable both the Prandtl and Karman mixing length results are shown. At 
the lower Mach number (M = 4) the differences between the various assumptions 
within a theory and between the theories is much less than at the higher Mach 
number. At M = 20 skin friction covers the gamut from increasing to 
decreasing with decreasing wall temperature. With one (the Van Driest) theory 
virtually all the possibilities can be obtained, depending upon the assumptions 
used in the calculations, including the case where the skin friction is prac- 
tically independent of wall temperature. This latter case is similar to 
results from the Spalding-Chi method though the numerical values are higher 
(even taking into account that the value of u> used in computing the Van 
Driest theory in figure 12 would have to be reduced in order to more closely 
correspond to the to of the data used by Spalding-Chi in formulating their 
result). Because of the sensitivity of this and similar theories to the vari- 
ous assumptions, there is more likelihood of definitive experiments being per- 
formed in the hypersonic speed range than at lower speeds. 


THE MAXIMUM HEAT TRANSFER TO SMOOTH FLAT PLATES AND CONES 


Up to this point the results do not allow an evaluation of the maximum 
heat transfer to be encountered on smooth flat plates and cones. This maximum 
is closely represented by the peak in heating which has been used to designate 
the virtual origin of the turbulent boundary layer in the previous section of 
this paper. Examples of the peaks shown by experimental heating distributions 
are contained in the upper portion of figure 1 and in figure 2(a) while in the 
upper portion of figure 4 such peaks are shown in both shear and heating 
distributions. 

Figure 13 has been prepared to show the values of peak heating (Ngt,p) 
obtained on flat plates on both low speed [j5 J and high-speed flow as a function 
of Reynolds number based on the distance from the leading edge to the peak 
heating location. The high-speed data is that used to compose figures 2, 3> 
and 4, plus unpublished data obtained by Sterrett in the Langley 20-inch hyper- 
sonic tunnel and data contained in (a 6] . The values of the peak in the heat- 
transfer coefficient show much the same variation with Reynolds number shown by 
the incompressible turbulent theory (Karman-Schoenherr). This suggests uti- 
lizing this theory to remove much of the Reynolds number dependence as is done 
in the lower portion of the data presentation in figure 13 for each speed 
range . In this form the Reynolds number dependence is indicated to be weak.; 
however, a definite dependence on Mach number is shown. 

A means for predicting these peak values would be desirable. Momentum 
matching, which was discussed in an early part of this paper, suggests itself. 
The result from momentum matching for incompressible flow using the Blasius 
equation for laminar skin friction, and the Karman-Schoenherr equation for tur- 
bulent skin friction is shown in the upper part of figure 13. The agreement 
with the data is surprisingly good considering that the theory assumes sudden 
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transition and the actual transition can he a gradual process as shown in 
figure 1. 

Application of this method to the compressible results shown in the lower 
part of figure 13 met with mixed results. Utilizing the Spalding-Chi formula- 
tion [20J with the Blasius equation modified to- compressible flow by the 
T-prime method gave good agreement with results but underpredicted the 

data by increasing amounts as the Mach number increased until at = 8 it 
was about 35 percent below the data. Some rough calculations with Li-Nagamatsu 
theory [25] for the turbulent flow gave values within 10 percent of the M2 = 8 

data. These calculations are not shown, as at present they require too much 

judgment for application (only the curves of figure 11 were available, limited 
to a Reynolds number of 10 ?, adiabatic wall, and average skin friction) . 

Another way of looking at peak data is to compare it with data downstream 

of the assumed virtual origin, which was shown previously to have a similar 

trend with Mach number. Figure lb presents these results in the form of the 
peak heating (or skin-friction) ratio divided by the heating (or skin-friction) 
ratio downstream of the virtual origin. The data is generally from sources 
cited earlier in this paper where the test conditions are given in some detail 
with the addition of Coles skin-friction data [2] taken under adiabatic wall 
and essentially room temperature supply conditions. To avoid a confusion of 
symbols, only average values are shown for the large number of peak results 
available from a given source. 

The flat-plate results for this ratio on figure lb while not conclusive 
appears to be independent of Mach number and the level agrees with that found 
for incompressible flow by the momentum matching procedure . For all the flat- 
plate data except the low-speed results the values of an d Ry are of 

the same magnitude. (The low-speed data were compared to results of momentum 
matching for the proper Reynolds number in figure 13 •) 

The same ratio from cone data is about 20 percent below the flat-plate 
results and appears also to be essentially independent of Mach number. (Note 
that the flight results are for Reynolds numbers about one order of magnitude 
higher than the other flat-plate or cone results shown.) Based on these results 
the prediction of the entire heat-transfer (or skin-friction) behavior on a flat 
plate or cone depends upon the ability to predict the Ng-fc/Nst,i ( or Cf/Cf^i) 
ratio . 


BOUNDARY-LAYER SURVEYS 


In the previous discussion, only gross effects of the turbulent boundary 
layers have been considered, i.e., skin friction and heat transfer. However, 
it is clear that a general understanding of the mechanisms at work within the 
boundary layer itself is needed for fully evaluating turbulent data and theories. 
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Some of the basic structure of the boundary layer can be obtained from conven- 
tional boundary-layer pressure and temperature surveys; however, because of the 
relative thinness of the laminar sublayer as compared to the overall boundary- 
layer thickness, small probes or thick boundary layers, or both, are necessary 
for probing in this region. 

Figure 15 is composed of recently obtained velocity and temperature pro- 
files in relatively thick boundary layers (2 to 4 in.) at Mach numbers of 6 and 
6.8 in air and 18.5 in helium. The ratios of wall to free-stream stagnation 
temperature are O.63, 0.50, and 1, respectively. This data was obtained on 
nozzle walls near the test section where the longitudinal pressure gradients 
were, for all practical purposes, zero. The Mach 6.0 data were obtained by Cary 
in the Langley 20-inch hypersonic tunnel, the Mach 6.8 data were obtained by 
Neal in the Langley 11- inch hypersonic tunnel, and the Mach 18.5 data were 
obtained by Watson in the Langley 22-inch helium tunnel. 7 The velocity profiles 
were obtained from local pitot pressure measurements and faired values of local 
stagnation temperature measurements® in conjunction with the conventional 
assumption of constant static pressure through the boundary layer. A recovery 
factor of unity was used for the temperature probes as time did not permit a 
comprehensive calibration of the probes; however, calibration of probes identi- 
cal to those used to gather the data, at Mach numbers of 6.0, 6.8, and 8.5> 
over a sizable range of unit Reynolds numbers gave recovery factors of unity. 

The recovery factor for these probes may or may not be unity at flow conditions 
encountered well within the boundary layer and especially near the wall. 

General examination of the velocity profiles in figure 15 show that they 
display the typical appearance of turbulent boundary layers, i.e., an outer 
region that can be represented quite well with a power law, an inner region, or 
laminar sublayer, and of course the intermediate, or buffer, region between the 
two. Closer examination of the profiles show that at a given Mach number the 
inverse of the exponent in the power law fit to the outer region increases with 
Rq, that there is a small inflection point near the edge of the boundary layer, 
that the laminar sublayer thickness (which is very small for the air data) 
increases with decreased Reynolds number, and that the values of Cf deducible 
from the velocity profiles could not be called conclusive. Walz [27] has criti- 
cized the validity of the Cf data of references [21J and [22] as well as 
Hill's older results [28] which were obtained from velocity profiles. Figure 16 
is composed of portions of representative velocity profiles from references [21J 
and [22] at a Mach number near 5* Included in these plots are results from 
reversed calculations using the values of Cf (and other pertinent quantities) 
given in the references and using the viscosity law given by equation (Al). In 

^Note on figure 15(c) that the boundary- layer thickness covers about half 
the radius of the nozzle, and thus significant lateral curvature effects are 
probably present. 

8ln general, the wall-distances at which the temperature measurements were 
made differed from those where the pressure measurements were made so faired 
values of the temperature data were used in computing the velocity profiles. 
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several instances, and especially at the lower values of T w /Tt,a> it is diffi- 
cult to determine how conclusive values of Cf could he deduced. The dashed 
curves, which were inversely computed from Spalding and Chi's [20] theoretical 
Cf values are shown here, and elsewhere on all the velocity profiles, except 
the one in helium, as reference lines. In many cases there is good agreement 
of experiment with the Spalding -Chi predictions; in others it is difficult to 
determine whether the disagreement is a fundamental one or a result of probe 
inaccuracies and interference effects near the wall. 

Returning to figure 15 for the temperature profiles, it is seen that 
Tf j Tf a, versus log y is roughly linear for the air data with some dependence 
on R e . The single profile taken in helium has a bowl-like shape, the lower 
value ^Tf |Tf^ OT ~ O.92) is probably close to the value that would be obtained 
at an insulated wall. There is some scatter in the M = 6.8 temperature pro- 
files which is believed due to changes with time in the previous history of the 
flow. During the gathering of this data it was observed that for a given y 
location, the measured values of Tf / Tf gradually increased with time (avail- 
able run time approximately 60 seconds). This probably resulted from sizable 
increases with time in the temperature of the nozzle wall upstream of the 
probing region. The test-section region, which is massive in structure and 
has relatively low heating rates, did not significantly increase in temperature 
during a run. Therefore, only temperature data taken near the start of the run 
(approximately 10 seconds) has been used since this data would correspond more 
nearly to a uniform wall temperature close to that measured in the probing 
region. The Mach 6.0 profiles, which were taken in a slow starting but long 
running time capability tunnel, and the Mach 18.5 profile, which was taken in 
a facility operating with a stagnation temperature near ambient temperature, 
were not found to be time dependent . 

Velocity and temperature profiles recently obtained by Adcock and Peterson 
on a hollow cylinder model in the Mach 6.0 Langley 20-inch hypersonic tunnel 
at T^yT.^,,0 values of 0.38, 0.4$), and O.89 are shown in figure 17- Qualita- 
tively, the data for the nonadiabatic wall temperature cases, agrees with the 
Mach 6.0 nozzle wall data previously presented in figure 15(a) except for the 
overshoot in T-t/T-t j00 . An overshoot in T-t/^t,oo for the adiabatic wall case 

is to be expected in order to conserve energy;^ however, this data shows a 
higher degree of overshoot for the nonadiabatic wall temperature cases than for 
the adiabatic one. Experimental difficulties with the small probes used to sur- 
survey the relatively thin boundary layers on the hollow cylinder possibly 
account for this behavior. By interpolating this hollow cylinder data for a 
T w /T t a, *= O.63 and extrapolating the nozzle wall data of figure 15(a) to lower 

values of Rg, it is seen that the resulting power in the power law fit to the 
outer u/U profiles do not agree. This result suggests that differences in 
past history of the flow are important even though the local conditions at the 
measuring station are the same. 

^From integration of the energy distribution across the boundary layer, 
conservation of energy is known to hold for Adcock and Peterson's adiabatic 
wall data. 
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In figures 18 and 19, stagnation temperature -velocity profiles in the cor- 
relation form (Tt - T w )/(T t , oo - T w ) versus u/U suggested by the Crocco 

energy equation are shown for the data previously presented against wall dis- 
tance as well as data from several additional sources [28], [29^ , [30] , [31] , 
and [32J . In this type of plot, Crocco 's equation for Np r = 1 gives a one- 
to-one correspondence between the two parameters and thus provides a convenient 
reference line to which the data can be compared. For cases where the y dis- 
tances at which the temperature data were taken differed from those at which 
the velocity data were computed, faired values of u/U have been used. Fig- 
ure 18 is composed of ^ = 6.0 data for several values of wall-to- stream stag- 
nation temperature ratio while figure 19 contains data at different Mach numbers 
for two wall temperature cases, the adiabatic wall case and the T w /T-t ,0° ^ 0*5 
case. From these figures, it is seen that significant deviations from the 
Crocco values do occur, depending upon the Mach number and wall temperature. 

For the constant Mach number data of figure 18, Walz' theory is seen to quali- 
tatively predict the trend of wall temperature, except for the overshoot in 
Tt/Tt^ near the outer edge of the boundary layer for Adcock and Peterson's 
data. As previously discussed, no overshoot in T-t/T-^oo would be expected for 
wall temperature significantly below adiabatic temperature , but an overshoot 
for the adiabatic case is expected. 

In the adiabatic wall-temperature profiles of figure 19, the flat-plate 
incompressible data of Reynolds, Kays, and Kline^O (29] are seen to be in excel- 
lent agreement with the Crocco equation. All of the compressible data show 
overshoots in Nothwang's flat-plate data J30J at M^, = 3-0 and 

Kistler's nozzle test section data (31] at Mach numbers of 1.7, 3-6, and 4.7 
lie mostly above the Crocco line while that of Lobb, Winkler, and Persh is bowed 
considerably below the Crocco line as are also Adcock and Peterson's Mach 6.0 
hollow cylinder data at intermediate values of u/U. Nothwang's and Kistler's 
data suggests a small downward trend with Mach number; however, there is a large 
unexplained difference between Kistler's Mach 4.7 data and the NOL Mach 4.9 data. 
Differences in past history of the flow or the fact that the NOL data were taken 
in a pressure gradient could be suspected for this behavior. However, as can be 
seen in figure 18 by visually noting the trend against wall-temperature ratio at 
the intermediate portion of the boundary layer, the Mach 6.0 hollow cylinder 
data (no pressure gradient) and the Mach 6.0 nozzle wall test section data (with 
previous pressure gradient history) seem to agree well with each other even 
though the velocity profiles were noted earlier to be significantly different. 
Also, the adiabatic wall hollow cylinder data at u/U values below about O.85 
and the Mach 4.9 NOL data agrees well with each other as well as Walz* theory, 
which is insensitive to Mach number. 

For the T v jT-t )00 0-5 profiles of figure 19, it is noted that the nozzle- 
wall data tends to be significantly lower than the flat plate or hollow cylinder 
data. Here, it thus appears that past history of the flow has a pronounced 
effect upon the profiles. Within the intermediate portion of the boundary layer, 

l^rhe values of Rq shown for this data were obtained using the Karman- 
Schoenherr equation and the given values of R x . 
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Danberg’s Mach G.k flat-plate data fp2j and Adcock and Peterson's Mach 6.0 
hollow cylinder data are in good agreement with each other as well as Walz’ 
theory. The NOL data [ 22 J , which covers Mach numbers from 5-1 to 8.2, shows 

no apparent effect of Mach number. The answer cannot be arrived at by the 
present superficial examination of the available data but will probably require 
a knowledge of the previous history of the boundary layer together with a 
sophisticated theory. 


CONCLUDING REMARKS 


Data from a large number of sources including flight and wind tunnel have 
been analyzed to provide information on the turbulent boundary layer at high 
speeds. Two of the obvious but important problems considered were the effect 
of Mach number and wall temperature on the heat transfer. In the supersonic 
speed range flight data had to be relied on to provide the trends at low wall- 
temperature ratios. The results were equivocal in that the data did not con- 
sistently support any one theory, though a trend of increasing heat transfer 
with decreasing wall temperature was indicated. For Mach numbers of about 5 
and greater, there was found to be little effect of wall temperature on the 
heat-transfer level as predicted by the Spalding-Chi method. There was, how- 
ever, a tendency for the Spalding-Chi method to underestimate the level of this 
heat transfer. There are indications that the Li-Nagamatsu theory, which takes 
into account the contribution of the density fluctuations to the shearing 
stress may account for the increased heating. The present experimental results 
do not allow a conclusive answer to this problem and definitive answers will 
probably have to await experiments at still higher Mach numbers than those pre- 
sented in this paper. 

The maximum heat transfer to flat plates and cones has also been examined. 
This maximum is closely represented by the peak in heating which has been used 
to designate the virtual origin of the boundary layer. For low-speed flow, 
momentum matching was found to give a good prediction of the experimental 
results on a flat plate. When properly ratioed to the results downstream of 
the virtual origin, the flat-plate peak heating ratios were apparently inde- 
pendent of Mach number at a value predicted for low- speed flow. The same ratio 
from cone experiments was lower than that given by the flat-plate experiments 
but also appeared to be essentially independent of Mach number. 

Lastly, the structure of the boundary layer has been considered. Attempts 
to correlate the data in the simple form suggested by Crocco’s integral of the 
energy-plus-momentum equation met with mixed results. For this type of plot, 
no definite trends with Mach number could be established, but a trend did 
appear to develop with wall -temperature variation which, in general, agreed 
with that predicted by Walz. Also, it was found that differences in past his- 
tory of the flow had a pronounced effect upon the profiles even though the 
local conditions at the measuring station were the same. 


15 


In regard to the determination of skin friction from the often used 
velocity profile method, it was noted that it is difficult to obtain conclusive 
results, mainly because of the few data points usually available in the thin 
laminar sublayer region as well as probe interference effects near the wall. 
Examination of some published profiles along with the deduced skin-friction 
coefficients revealed that several of the skin-friction values could be sig- 
nificantly in error due to these difficulties. 



APPENDIX A 


BASIC EQUATIONS USED IN THE PAPER 


Certain equations used in reducing and evaluating the data in this paper 
are included here for convenience and clarity. 


Viscosity of air: 

Keyes ' proposed three constant empirical formula [33] was used 

= 0.0232 T 3 / 2 x 10~ 6 lb-sec 
T + 220 x 10-(9/ T ) ft 2 

Karman-Schoenherr incompressible skin friction formulas: 

0.242 




- lo §io ( C F, i R x) “ log lo( 2 R e) 


= 


0.242C f> . 


0.242 + 0.8686^5^7 
The Karman form of Reynolds analogy factor [7] : 


-1 


2N 

I 

c f,l 


= S 1 + (N Pr - l) + log, 


5N pr + 1 


(Al) 


(A2) 


(A3) 


(A4) 


The application of equation (A4) to equation (A3) to obtain Ng^ ^ is straight- 
forward but the application to compressible flow is not. In applying the 
Karman factor to compressible flow the particular theory was transformed to the 
incompressible plane. This was possible with the T-prime and Spalding-Chi 
methods. Then the Karman factor corresponding to this incompressible (by 
transformation) value of Cf was used to change the compressible value of Cf 

to the compressible value of Stanton number. This was checked by plotting the 
Deissler-Loeffler [8j values for Reynolds analogy factor against the values of 
Cf transformed to incompressible flow by the method of Spalding-Chi (i.e., 
using the Spalding-Chi values of F c and the Deissler-Loeffler value of Cf 
to obtain F c Cf where F Q Cf is taken to be Cf^p). In this form the Deissler- 

Loeffler values of Reynolds analogy factor correlated remarkably well and were 
within 2 to 3 percent of the value given by Karman. In applying the Karman 
factor in this paper Np r was taken as O.725. 
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APPENDIX B 


TRANSFORMATION OF THE LOCAL SKIN FRICTION ON A 
TRUNCATED CONE TO THAT ON A POINTED CONE 


If the local skin friction coefficient may be represented by the power law 
equation Cf « R _1 / n then the transformation of distance on an axisymmetric 
body to that on a corresponding two-dimensional shape is given by Mangier [ 31 +] 


x 



(Bl) 


y = 



(B2) 


The bar designates the quantities in two-dimensional flow; r Q is the local 

radius of the axisymmetric shapes; x is measured along a geodesic on the body; 
y is on a normal from the body surface. With n = 2, corresponding to laminar 
flow, the familiar exponent of 2 is obtained in equation (Bl). 

From equations (Bl) and (B2) and the momentum equation one obtains the 
ratio of the local skin friction on a three-dimensional body to that in two- 
dimensional flow as 


1 1 


C f Ry 1 / 11 

/ \ n ~ 1 
cro) 

X 


Cf Rv 1 / 11 

\ L I 

r x /cv 

J 0 V d , 

n 

in 

) dx 


(B3) 


For the specific case of the ratio of the local skin friction on a truncated 
cone to that on a flat plate the result is 
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or the ratio of the local skin friction on a truncated cone to that on a 
pointed cone is 



Equation (B5) with n = k was used to correct the turbulent cone data to 
correspond to that on a pointed cone with turbulent flow from the apex. The 
sketch below defines the Reynolds numbers shown in equations (B4) and (B5). 
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Figure 5-- Heat-transfer coefficients on a flat plate in a high-speed flow as a 
function of Reynolds number measured from peak heating. Data from Langley 
l8-inch variable-density wind tunnel; d = 0.002 in. ( 0.05 mm). 
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Figure 8.- Heat-transfer-coefficient ratio as a function of wall temperature 
ratio from flight and wind-tunnel experiments in supersonic flow and as 
given by two prediction methods. 
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Figure 10.- Heat- transfer coefficient ratio as a function of wall temperature (enthalpy) rati 

in hypersonic flow. 
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Figure 11.- Local and average skin-friction ratios as a function of Mach num- 
ber as predicted by various theories j Reynolds number of 10?. 
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Figure 12.- Various predictions of the effect of wall temperature on local skin 

friction; Reynolds number of 10?. 
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skin-friction) ratio downstream of the virtual origin as a function of 
Mach number. Shaded symbols - skin-friction dataj open and solid symbols - 
heat-transfer data. 
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(a) Moo = 6; profiles on the test- section wall of the Langley 20-inch 
hypersonic tunnel (two-dimensional contoured nozzle). 

Figure 15*- Velocity and temperature profiles obtained in thick turbulent 
boundary layers growing on nozzle walls. 
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Figure 17*- Velocity and temperature profiles at = 6 obtained on a hollow 
cylinder model in the Langley 20-inch hypersonic tunnel. 







Figure 18.- Effect of vail temperature on the temperature-velocity profiles at 
= 6 in the Langley 20-inch hypersonic wind tunnel. 



N1 


x 10 


-4 


Model 


Source 


• 0 

■ 0 

O 1.7 
□ 3.6 
O 4.7 

//////3.0 


.17 \ 

.45 

4 .00 ^ 
3.30 > 
2 . 88 ^ 
1.34 


1.03 535 Flat plate 


540 

525 


Noz. test sec. 
Flat plate 


Reynolds, Kays, 
and Kline 

Kistler 

Nothwang 



® 6.0 
^ 6.4 
■ 6.8 


1.15 0.49 865 Hollow cyl. Adcock and Peterson 


.60 

1.30 


.52 

.50 


990 

1110 


Flat plate Danberg 
Noz. test sec. Langley 11 in. 



NASA 


Figure 19.- Temperature-velocity profiles for two wall temperature ratios and 
various values of Mach number at the edge of the boundary layer. 
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